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ABSTRACT 

We propose a new method of estimation of the black hole masses in AGN based on the 
normalized excess variance, Cj^xs- We derive a relation between ct^xsj the length of the 
observation, T, the light curve bin size, Ai, and the black hole mass, assuming that 
(i) the power spectrum above the high frequency break, fbf , has a slope of —2, (ii) the 
high frequency break scales with black hole mass, (iii) the power spectrum amplitude 
(in frequency x power space) is universal and (iv) cr^^g is calculated from observations 
of length T < l/i^hf- Values of black hole masses in AGN obtained with this method 
are consistent with estimates based on other techniques such as reverberation mapping 
or the MeH-stellar velocity dispersion relation. The method is formally equivalent to 
methods based on power spectrum scaling with mass but the use of a^^^ has the big 
advantage of being applicable to relatively low quality data. 
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1 INTRODUCTION 



The X-ray emission of active galactic nuclei (AGN) displays 
variations over a wide range of time scales. The first con- 
vincing demonstration of this phenomenon came with the 
EXOSAT 'long looks' (Lawrence et al. 1987; McHardy & Cz- 
erny 1987). The data showed no characteristic time-scales, 
and the power spectral density function (PSD) showed a 
'red-noise' shape. At the same time, it was also noticed 
that more luminous sources show slower variations. Various 
methods have been used in the past in order to determine 
the X-ray variability amplitude of AGN: determination of 
the 'two-folding' time-scale (i.e. the time-scale for the emit- 
ted flux to change by a factor of two), calculation of the 
PSD amplitude at a given frequency, and estimation of the 
so called 'normalized excess variance' (cTnxsi i-e. the variance 
of a light curve normalized by its mean squared after correct- 
ing for the experimental noise). In all cases, these quantities 
appear to anti-correlate with the source luminosity (Barr 
& Mushotzky 1986; Lawrence & Papadakis 1993; Green, 
McHardy & Lehto 1993; Nandra et al. 1997; Turner et al. 
1999; Leighly 1999; Markowitz & Edelson 2001). One pos- 
sible explanation for the observed anti-correlation between 
the X-ray variability amplitude and X-ray luminosity, Lx, 
is that more luminous AGN have larger black hole masses 
(Mbh) as well. In this case, as Mbh increases, so does the 
size of the X-ray source, and a change of the source lumi- 
nosity by a constant fraction takes relatively longer; equiv- 
alently, variability amplitude measured at a fixed timescale 



should decrease with mass. The significant progress in the 
measurement of the black hole (BH) mass in the centers 
of nearby galaxies has made it possible to actually directly 
test the dependence of the X-ray variability amplitude on 
Mbh in AGN. The data show that cr^xs measured at a cer- 
tain timescale is indeed anti-correlated with Mbh (Lu & Yu 
2001; Bian & Zao 2003; Papadakis 2004). 

If there is an intrinsic correlation between X-ray vari- 
ability and Mbh , then X-ray variability measurements could 
be used in order to measure the central black hole mass 
in these objects. This possibility has already been stud- 
ied in the past. Hayashida et al. (1998) and Czerny et al. 
(2001) have used the PSD normalized by the square of the 
mean flux as a measure of the X-ray variability of a source. 
By calculating the ratio of the frequencies at which the 
PSD X frequency has a certain value in AGN and Cyg X- 
1, they were able to estimate the Mbh in AGN. Recently, 
long, well-sampled RXTE light curves have been used in or- 
der to accurately estimate the X-ray PSD of AGN. One of 
the main results from these studies has been the unambigu- 
ous detection of a characteristic 'break frequency', i^bf, at 
which the power spectrum changes its slope from a value of 
—2 to —1 above and below Vhf, respectively (Uttley et al. 
2002; Markowitz et al. 2003). This break frequency is anal- 
ogous to the high frequency break in galactic sources and 
should not be confused with the low frequency break, where 
the spectrum changes slope from —1 above to below. This 
high frequency break appears to correlate well with Mbh, in 
the sense that oc Mbh (Markowitz et al. 2003). Thus 
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another way to estimate Mbh is through the determination 
of fbf in the X-ray PSD of a source. However, this method 
requires rather long, high quality data. 

We propose a new method to estimate Mbh in AGN, 
which may be used even for relatively low quality and/or 
short data coverage. The method is based on the use of cr^xs 
as a measure of their X-ray variability amplitude, and takes 
Cyg X-1 as a reference point. Using recent, archival RXTE 
and ASCA light curves of 21 AGN, we estimate their Mbh, 
and we compare our results with Mbh estimates which are 
based on other methods. We find that the accuracy of the 
Mbh estimates using our new method is roughly similar to 
the accuracy of the estimates with well established methods 
like reverberation mapping. Furthermore, the a^xa estima- 
tion of an X-ray light curve is rather easy and straightfor- 
ward. Gonsequently, as long as typical AGN observations 
performed by present day satellites (like XMM-Newton and 
CHANDRA) detect significant variations, they can be used 
to estimate (Jnxs and hence Mbh, using the method which 
we present in this work. It is possible then that the use of 
our method can provide Mbh estimates for a large number 
of AGN in the near future. 



2 THE RELATION BETWEEN BH MASS AND 
EXCESS VARIANCE 

The scaling of the PSD with mass is most conveniently dis- 
cussed after multiplying the normalized PSD by the fre- 
quency ly. We assume here that the break frequency, vm, 
scales with the mass: fbf oc Mbh~^, but the amplitude of 
the normalized PSD x v part at the frequency break does not 
depend on mass. We specifically define this normalization at 
Vhi and call it the PSD amplitude. 



PSD^mp = -P(l'bf) X i/bf. 



(1) 



where P{iy) is the power spectrum normalized to the mean 
squared. Recently, Papadakis (2004) has showed that this 
quantity is indeed similar in a few well studied Seyfert 1 
galaxies. 

Let us suppose that an AGN is observed for a period 
T < 1/vhi, and that the resulting light curve is binned using 
a bin size of At. The excess variance of the light curve is 
equal to the integral of the PSD, 



l/2At 



P{u)du 



1/T 



Let us also suppose that P(«^) 



(2) 



where 



A =PS-Damp/t'bf and PSDiimp is constant, and that vm 
— _B/Mbh, where B is another constant. Under these as- 
sumptions, 

B 



'^nxs — PS Damp ^ 



-(T - 2 At). 



(3) 



The equation above can be re-written as follows, 
Mbh =C{T -2 At) /aL„ 



(4) 

where C = PSD^mpB. Eq. can be used to estimate Mbh 
provided the a^^g of a source has been estimated, and C is 
known. In order to determine C, we used the results from 
the PSD analysis of a Galactic black hole candidate, namely 
Cyg X-1. 



The X-ray PSD of this source has been extensively 
studied the last twenty years or so (Kemp et al. 1978; 
Gierlinski et al. 1999; Nowak et al. 1999; Ghurazov et al. 
2001; Pottschmidt et al. 2003). In its low/hard state, the 
PSD of the source can be roughly represented as a power 
law with a high frequency slope of —2 , which breaks to a 
slope of —1 below 2 — 3 Hz, and then to a slope of zero be- 
low ~ 0.2 - 0.02 Hz (Nowak et al. 1999; Pottschmidt 2003). 
Assuming that Eq. I^J holds for Gyg X-1 as well, and adopt- 
ing the value of 10 Mq for its black hole mass (Nowak et 
al. 1999; Gierlinski et al. 1999 and references therein), we 
determined the value of C as follows. 

We used the resuhs from the PSD analysis of 130 RXTE 
observations of Gyg X-1 reported by Pottschmidt et al. 2003. 
They calculated the PSD for each observation and fitted 
them using Lorentzian functions. The RXTE observations 
cover a long period of time (from December 1997 to July 
2001) in which Gyg X-1 was mainly in its hard state. We 
chose 68 observations, carefully avoiding data from soft or 
failed transition states. We calculated the excess variance 
using the relation ''"nxs^/j^^^ P{i/)du, where the lower fre- 
quency limit was chosen in such a way that it will always be 
higher than //bf in this source. Using Eq. (4), with Mbh = 10 
Mq, At = 0, and T = 1/(10 Hz)= 0.1 sec, we find that C 
varies between ~ 0.7 and ^ 1.3, with the arithmetic mean 
being C — 0.96 ± 0.02. We accept this value for the normal- 
ization constant C in Eq. I^J. 



3 NEW Mbh ESTIMATES FOR 21 AGN 

In order to investigate the accuracy of our method, we 
selected for analysis several AGN with known black hole 
masses and recently obtained RXTE or ASCA light curves. 
First, we estimated the (Tnxs of each light curve using the 
relation (Nandra et al., 1997), 



A'"data 



E 

i=l 



(5) 



where A^data is the number of data points, and x, (Ji are the 
unweighted, arithmetic mean and error of the Xi, respec- 
tively. Then, using Eq. (|1J we estimated Mbh, and compared 
our Mbh estimates with the existing values from literature 
based on the reverberation mapping technique or the mea- 
surement of stellar velocity dispersion (Section l3.2l Table 1). 



3.1 Remarks about preparation of the data 

RXTE has regularly observed a large number of AGN over 
the last few years, providing long light curves with good 
signal to noise (S/N). For almost all the objects listed in 
Table 1 we have retrieved 2 — 9 keV PGA background sub- 
tracted light curves from the standard data products ('Std- 
Prods') archive of RXTE (detailed information about the 
data in this archive is given at the RXTE mission web 
site). There are no archival RXTE data for 3 AGN with 
previously estimated Mbh (namely PG121H-143, NGG4395 
and NGG4593). For these sources we retrieved 2 — 10 keV, 
ASCA light curves from the TARTARUS database. We used 
background subtracted GIS and SIS light curves which we 
combined in all cases, except from the first observation of 
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Table 1. Previous and new Mbh estimates for 21 AGN 



Name 


Mbh,l 




Proposal 




(xIO'^Mq) 


(xlO'^M©) 


number 


3C120 


6.3 (1,2) 


14.5 


R /P30404('3") 


3C390.3 


35.5 (1) 


26.8 


R/P10340(7) 


Aknl20 


18.6 (1) 


5.6 


R/P30232(l) 


IC4329A 


6 fll 


12.3 


R/P40153(2) 








R/P20315(2) 








R/P10313(4) 








R/P50706(l) 


IC5063 (32) 


5.5 (2) 


(0 19) 

3.8 


R/P10337(l) 


Mrk34S ('S2"l 


l.V (^) 


7.5 


R /PI 0326('3") 


Mrk509 


7 5 CD 


7.0 


R/P10311(8) 


Mrk766 (*) 


0.35 (2) 


(0.0066) 


R/P60135(23) 






0.132 




NGC3227 


4:1 (1 31 


2.2 


R/P40151(6) 








R/P10292(3) 


NGC3516 


1 7 (3) 


1.7 


R/P50159(l) 








R/P20316(4) 








R/P30224(2) 


NGC3783 


1 (1) 


0.83 


R/P10297(2,l) 








R/P30227(l) 


NGC4051 (*) 


0.05 (5) 


(0.0014) 


R/P50153(128) 






0.028 




NGC4151 


1 4 (1) 


1.9 


R /P501 55('2") 








R/P00022(3) 








R/P00024(4) 


NGC4395 


0.007(6) 


0.022 


A/76006000(2) 


NGC4593 


0.66 (3) 


0.60 


A/71024000(3) 








A/75023010(l) 


NGC5506 (*) 


8.8 (4) 


(0.32) 
6.4 


R/P20318(10) 


NGC5548 


9.2 (1,7) 


12.3 


R/P30220(2) 








R/P10297(3) 








R/P30218(3) 


NGC7469 


0.7 (1) 


1.3 


R/P10315(20) 


PG0052+251 


26.1(1) 


10.0 


R/P40157(l) 








R/P20338(l) 


PG0804+761 


7.6 (1) 


5.6 


R/P40157(3) 


PG1211+143 {*) 


3.2 (1) 


(0.40) 


A/70025000(l) 






8.1 





Col. (1) lists the object name and type: (S2) denote Seyfert 
2 objects, (*) - Narrow Line Seyfert 1. All other objects are 
classified as Seyfert 1 or quasars. Col. (2) lists the Mbh l values 
taken from literature. The numbers in parentheses in Col. (2) 
correspond to the following references: (1) Kaspi et al. (2000), 
(2) Woo & Urry (2002), (3) Onken et al. (2003), (4) Papadakis 
(2004), (5) Shemmer et al. (2003), (6) Filippenko & Ho (2003), 
(7) Peterson & Wandel (2000). Col. (3) lists the Mbh,^^ values 
obtained in this work. For NLSls, the values in parentheses are 
the Mbh values found using Eq. (3), while the values listed 
below are multiplied by a factor of 20 (see text for details). 
Col. (4) lists the observation details. The first letter refers to the 
satellite (R - RXTE, A - ASCA), the number that follows refers 
to proposal number, and we list the number of light curve parts, 
A'^p, which we used in order to determine cr-^-^^, in parentheses. 

NGC4593. In this case we used the GIS data only, as there 
were many missing points in the SIS light curve. 

The character of the obtained light curves depended 
mostly on the type of the monitoring. A number of the light 
curves came from sparse extensive monitoring. The source 
was observed typically for 10 ksec once per day, for several 



days. For example, 3C390.3 was observed usually once per 
day, for 2 months, with a single observation lasting about 
3 ksec. Such light curves were binned to 5400 sec (roughly 
the orbital period of the satellite). If the resulting light curve 
consisted predominantly of single points separated by a day, 
then the remaining few consecutive data points were aver- 
aged to create a uniformly covered light curve. 

Other light curves came from intensive monitoring. 
The source was observed continuously for 40 ksec or more, 
with gaps only caused by earth occultation. For example, 
NGC3516 was observed for over one day. The observa- 
tional campaign was sometimes repeated or supplemented 
by sparse extensive monitoring. NGC7469 was almost con- 
tinuously observed for about 1 month: 10 to 20 ksec expo- 
sures were made almost every day. Such light curves were 
also binned to 5400 s. 

Several sources were occasionally observed for 10 — 
40 ksec. The bin size used for these light curves typically 
varied from 512 s to 5400 s, depending on the length of the 
data and the S/N ratio. 

The light curve of PG0052-I-251 was binned to 10800 s 
in order to increase the S/N ratio. The light curve from 
continuous monitoring of NGC4593 was binned to 2048 s 
and divided into separate sets of 30 ksec duration in order 
to estimate (Jnxs in more than 2 sub-parts. Even shorter bin 
size was requested for NLSl sources, as discussed below. 

The long light curves were divided into separate shorter 
strings of length T, as required by the condition T < 1/vhi 
(see Section 2). The value of Vhf for each source was es- 
timated from the previous Mbh estimates, using the 
Mbh relation of Markowitz et al. (2003). The PSD break fre- 
quencies of 3 NLSl objects (NGC4051, NGC4253/Mrk766, 
and NGC5506) are high (see McHardy et al. 2003; Vaughan 
& Fabian 2003 and Uttley et al. 2002). Consequently, their 
light curves were divided into parts with T =1.25, 2, and 
10 ksec, respectively. In order to have as many points as 
possible in them, we used a bin size of 16 s for the NGC4051 
and Mrk766 light curves, and a bin size of 128 s in the case 
of NGC5506. In some cases, even if T < l/t'bf, long, well 
sampled light curves were also divided into shorter parts in 
order to increase the number of ct^xs estimates. 

The cr^xs values were calculated for each short light 
curve or each substring of a long light curve separately. The 
number of light curve sub-parts (A'p) used for the Mbh es- 
timation of each source is listed in column (4) of Table 1. 
In the case when A'^p > 1, we used the average cTnxs value to 
derive the black hole mass from Eq. Q. 

3.2 Results 

The values of the black hole masses obtained from Eq. 
are given in Table 1. 

Column (1) in this table lists the names of AGN. 
Four objects are classified as 'Narrow Line Seyfert 1' 
(NLSl), while the others are classical broad line Seyfert 1 
(Sl)/quasars (QSO), and Seyfert 2 (S2) objects. Column 

(2) gives the Mbh estimates taken mostly from Kaspi et 
al. (2000), and they correspond to the mean value of their 
'mean' and 'rms' values. In cases where we found more than 
one Mbh estimate for an object, the value listed in Table 1 
corresponds to the mean of the respective estimates. Column 

(3) gives our mass determination based on X-ray variability. 
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Log Mn„, L (M^) 

Figure 1. Black hole mass estimates using the method presented 
in this work (Mgjj ^2 ) plotted versus the literature estimates 
(Mbh l)i for the objects listed in Table 1. Solid squares and solid 
triangles show the classical Seyfert 1/Seyfert 2/quasars and the 
NLSl estimates, respectively. Open triangles and the open square 
show the NLSl and IC5063 estimates after they have been multi- 
plied by a factor of 20 (see text for details). The solid line shows 
the Mgjj ^2 = Mbh,l relation. 



Fig. 1 shows a comparison of the two mass de- 
terminations. Generally, the agreement between Mgjj ^2 
and Mbh.l is quite good for the classical Seyferts (filled 
squares), but not in the case of NLSls (filled triangles). 
Our Mbh estimates for the latter are significantly smaller 
than Mbh,l. Excluding NLSls, we find that on average, 
(MBH,a2„/MBH,L) = 2.44±1.16. Two Other sources with the 
largest differences between Mbh ^-^d Mbh.l are IC5063 
and IC4329A; excluding also those sources from our sample 
we find that on average, (Mbh.ctS /Mbh.l) = 1.41 ± 0.31. 

The results for NLSl galaxies can be reconciled with 
the previous mass determination if we multiply our Mbh 
NLSl estimates by a factor of 20 (open triangles in Fig. 1). 
Interestingly, when we multiply our IC5063 estimate by the 
same factor 20 as for other NLSl, the new value becomes 
very similar to Mbh,l for that object. This source is classi- 
fied as S2 with broad polarized Balmer lines (Veron-Cetty 
& Veron, 2001). Colina et al. (1991) found strong evidence 
for the presence of a hidden, luminous source in the nucleus 
of IC5063. Perhaps then, this S2 galaxy may host a NLSl 
rather than a SI nucleus, like NGC5506 (Nagar et al. 2002). 
In this light we suggest that IC5063 is a hidden NLSl. 

The source of the problems with IC4329A is more diffi- 
cult to understand. If we use the stellar velocity dispersion 
measurement of 225 km s~^ (Oliva et al., 1999) for this 
object, together with the MBH-stellar velocity dispersion re- 
lation of Tremaine et al. (2002), we obtain Mbh~ 2.2 x 10** 
M0, entirely consistent with our estimate for this source. 
Marziani et al. (1992) found that this source has large ex- 
tinction and they gave E{B — V) ~ 0.8. After correcting the 
optical luminosity, vLi,, at 5100 A for the large extinction, 
Nikolajuk (2004) finds Mbh= 3.9x10® M©, using the accre- 
tion disk method. This value is again consistent with the 



Mbh estimate listed in Table 1. Therefore, the value given 
by Kaspi et al. (2000) may be too small. 

When we exclude IC4329A from our sample but we con- 
sider SI mass estimates, together with the scaled NLSl and 
IC5063 estimates, then (MBH,a2„/MBH,L> = 1.30 ± 0.25. 
This result suggests that the new method we propose for 
the Mbh estimation in AGN provides estimates which are 
consistent with the estimates based on the reverberation 
mapping technique or the MBH-stellar velocity dispersion 
relation. 

4 DISCUSSION 

In this work we propose a new method for the estimation of 
black hole mass in AGN. The method is based on the use of 
the 'normalized excess variance', cTnxsi ^ estimated from X- 
ray light (Eq. |KJ. We find a simple expression between (Jnxs 
and Mbh (Eq. 4), if we assume that the power spectrum in 
AGN has a 'universal' form: (1) there is a break frequency 
which scales with Mbh, (2) the PSD slope above the break 
frequency is —2, and (3) the PSD amplitude below the break 
frequency is roughly constant in all objects, independent of 
Mbh (see Section 2). Note that the method is applicable to 
2—10 keV light curves, as the PSD may be different in other 
energy bands (e.g. Nandra & Papadakis 2001; Vaughan & 
Fabian 2003; Vaughan et al. 2003). 

Using our method on archival RXTE, and ASCA, 2-10 
keV light curves, we estimated Mbh in 21 AGN, with pre- 
viously known Mbh estimates. The comparison between the 
new and past Mbh estimates indicates that the method 
works very well for 'classical' AGN, i.e. SI, S2 galaxies and 
quasars. The Mbh values that we find with the proposed 
X-ray variability method are comparable to the values de- 
termined with other methods. However, there seems to be 
a significant disagreement for NLSl objects. Their mass 
estimates, determined from the X-ray variability method, 
are significantly smaller than the Mbh estimates which are 
based on other methods. The basic issue, that there is a 
shift in the X-ray variability amplitude of NLSl and 'classi- 
cal' AGN of the same mass, was already noticed by Czerny 
et al. (2001). Furthermore, Papadakis (2004) and McHardy 
et al. (2003) have showed that the break frequency Vhi of 
NLSl is probably higher than the respective frequency in 
the power spectra of 'classical' AGN with the same black 
hole mass. If B^^^^ ~ 20 x B'^'^'^^ (Papadakis 2004), where 
B is the constant defined in Section (2) , and PSDamp is the 
same for both classes of objects, then since C =PSDampB, 
we should expect that C^lsi ^ 20 x C"^'''"^. If we adopt this 
view, then we obtain an excellent agreement between the 
Mbh estimates which are based on the present X-ray vari- 
ability method and the previous estimates for all objects, 
both 'classical' AGN and NLSls. 

The method we present here is based on the assump- 
tion of a 'universal' power spectrum shape in AGN. Con- 
sequently, it is formally equivalent to other methods which 
directly use the scaling of a certain power spectrum char- 
acteristic (like the break frequency, or a fixed PSD value) 
with Mbh (e.g. Hayashida et al. 1998; Czerny et al. 2001; 
Markowitz et al. 2003). However, as our method requires 
the estimation of o"nxs only, it should be applicable for many 
more cases than methods which require the estimation of 
the PSD, as the latter is much more complicated. 
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Present day satellites like XMM-Newton and CHAN- 
DRA will soon provide a large number of AGN light curves 
with a typical exposure time of 30 — 50 ksec. Due to the low 
internal background of their detectors, most of these light 
curves should be of high signal to noise in the 2 — 10 keV 
band. Furthermore, because of their long orbital periods, 
light curves (as long as ~ 1.5 days) will be evenly sampled. 
Provided that significant variations are observed, they could 
be used to estimate (Jnxsy a-nd hence Mbh through Eq. (4). 

The choice of the bin size. At, is not crucial and should 
depend mainly on the average count rate. It needs to be large 
enough in order to minimize the effects of the experimental 
Poisson noise, but at the same time the resulting number of 
the points in the light curve, A'data, should also be large (at 
least larger than 10). 

The choice of the length of a single data string, T, must 
satisfy the requirement 1/T > ^bf- If i^hi for a given object 
is not known in advance, this condition should be checked 
a posteriori. For NLSl, the length should be ~ 20 times 
smaller than the appropriate T value for classical AGN, ac- 
cording to the results of the present work. 

Apart from observational errors, the determination of 
the CTnxs is biased by intrinsic statistical error connected 
with the specific shape of the power spectrum (Vaughan et 
al. 2003). Averaging over several independent variance mea- 
surements performed at various epochs reduces the error. If 
we have a single long light curve the error reduction is less 
efficient but can be assessed through simulations. We per- 
formed such simulations for NGC 5506 light curve (6 times 
longer than l/vhf) with many gaps, and with typical dura- 
tion of a single string of 10* s. Systematic 'leak' of the power 
from long to shorter timescales appeared to be negligible, 
and the intrinsic statistical scatter was A log o-nxs=lo i9' 
90 per cent confidence level. 

Finally, we determined the constant C in Eq. ||1J as- 
suming that Mbh = 10 M© for Cyg X-1. However, if a 
higher value of ~ 20 Mq is more appropriate for this ob- 
ject (e.g. Ziolkowski 2004), then C and the Mbh estimates 
listed in Table 1 should be increased by a factor of ~ 2. 
Such a change would, however, actually introduce a discrep- 
ancy between the masses estimated from the present X-ray 
variability and from other methods. Although the mass de- 
termination by the reverberation method may also contain 
a systematical error of the order of a factor 2 — 3 (see Kro- 
lik 2001), this possibility underlines the importance of the 
accurate determination of the properties of Cyg X-1 as this 
object is routinely used as a reference source. 



5 CONCLUSIONS 

• We propose a new method based on measurements of 
a^xs of X-ray variability of AGN to Mbh estimates. The 
estimates of black hole masses in AGN obtained from our 
method (Mgjj ^2 ) are consistent with the estimates based 
on other techniques such as the reverberation mapping tech- 
nique or the MBH-stellar velocity dispersion relation. 

• The constant C in Eq. @ is equal to 0.96 for 'classical' 
AGN and 19.2 (i.e. 20 times larger) for NLSl 

• We suggest that IC5063, classified as S2 galaxy, is in 
fact a hidden NLSl. 
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